Figure 1: The photon energy spectrum of B — > ~{X S process calculated by 
the parton model. 



Figure 2: Peterson's fragmentation function and our refined form for de- 
scribing decay processes. The solid line denotes our refined function and the 
dashed line Peterson's original function. 



Figure 3: (a) The photon energy spectrum of B — > ~{X S process with Peter- 
son's function (dashed line) and with our refined function (solid line), (b) 
The electron energy spectrum of B — > evX u process with Peterson's function 
(dashed line) and with our refined function (solid line). 



Table 1: The second moments of line shape a 2 , average kinetic energies of 
6-quark, effective masses of light degrees of freedom and 6-quark masses 
calculated with our refined function varying the parameter olq. 
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Abstract 



We test the parton model approach to inclusive B ~fX s and B evX u 
processes by calculating a few moments of the distribution function. The 
1/niQ expansion of 6-quark distribution function is discussed. We show that 
Peterson's fragmentation function is useful as the distribution function of b- 
quark but requires some improvements. By extension of Peterson's arguments, 
we obtain the refined form for the distribution function and with the refined 
function, we obtain the consistent results with the HQET. 
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I. INTRODUCTION 



The parton model approach has been established as a phenomenological model to de- 
scribe inclusive semileptonic B decays [1,2] and found to give good agreement with experi- 
ments for electron energy spectrum at all energies. The parton model is motivated by the 
deep inelastic scattering process (DIS) and it is well known that inclusive semileptonic de- 
cays of heavy flavours are intimately related to DIS via channel crossing. Jin and Paschos 
gave a theoretical foundation on this model recently [3]. The light-cone dominance of B 
decay processes allows one to use the methods of DIS and express the commutator of two 
currents as bilocal operator of quark fields. Its Fourier transform is related to the distribu- 
tion function of 6-quark inside B meson. Nonperturbative effects of QCD are included in 
the matrix elements of the bilocal operator and also in the distribution function. 

Over the past few years, there has been much progress in the study of inclusive semilep- 
tonic decays of B mesons with the help of heavy quark effective theory (HQET) [4-6]. One 
can get a QCD-based expansion in powers of l/rriQ by performing an operater product 
expansion on the hadronic tensor. With this technique, the systematic treatments of the 
inclusive spectra is achieved in principle. Therefore it is important to reanalyze phenomeno- 
logical models in the framework of the HQET. 

The ACCMM model [7], which is the most popular model for the inclusive semileptonic 
decays of heavy mesons, has been reanalyzed in the framework of the HQET recently by 
several authors [8,9]. They showed that ACCMM function for certain values of parameter 
can be viewed as good approximation for leading twist nonperturbative contribution to 
differential distribution curves. 

In this paper we analyze the parton model with the help of the HQET. Since the parton 
model is constructed in the analogy of DIS, we can study the model in the context of 
the HQET more easily than the ACCMM model. With properly defined mj, the parton 
model is represented by 1/m^ expansion, and the distribution function of 6-quark inside a 
B meson used in this model is directly related to the genuine QCD distribution function 
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which is not obtained systematically. Firstly we use Peterson's fragmentation function as 
the distribution function of 6-quark following the original works of ref. [1,2]. QCD cannot 
determine the genuine distribution function completely, but only parametrize its moments. 
Thus we represent the function as infinite sum of singular functions in power of l/rriQ with 
the moments. By calculating the moments of the QCD distribution function, we test the 
validity of the model. We find some discrepancies in the values of the parameters fj 2 v and A 
obtained by the parton model from those of the HQET while the parton model gives well 
agreed decay spectra with those of the HQET. Hence an improvement of the functional form 
is required. We refine the functional form by extending Peterson's approach to the decay 
process and obtain the consistent result with the HQET. 

This paper is organized as follows: in section 2 we formulate l/rriQ expansion of the 
distribution function. In section 3, we derive the differential decay rate of inclusive B — > 
~{X S decay in the parton model and show that the 6-quark distribution function of the 
parton model is corresponding to the genuine distribution function of QCD. The inclusive 
semileptonic b — > u decays are also studied. We test the distribution function by calculating 
a few moments in section 4. In section 5, we improve the functional form of Peterson's 
function. Our conclusions are summarized in section 6. 

II. HEAVY QUARK EXPANSION OF THE PARTON MODEL 

The parton model approach suggested in ref. [1,2] has been established on the bases of 
following two postulations: 

i) This model pictures the mesonic decay as the decay of the partons with the incoherent 
assumption in analogy to deep inelastic scattering process. The probability of finding a b- 
quark in a B meson carrying a fraction x of the meson momentum in the infinite momentum 
frame is given by the distribution function f(x). The decay probability is determined by 
the momentum of 6-quark inside the B meson and has the value of f(x)dx. Then we write 
the Lorentz invariant decay width as follows: 
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E B dT(B X q ev) = J dx f(x) E h dT(b qev) , (1) 
with the relation p b = xp B . 

ii) Peterson's fragmentation function [10] is used as the distribution function f(x) of 6-quark 
because we expect that the functional forms of the distribution function and fragmentation 
function are similar for heavy flavours. This is the fragmentation function usually used in 
Lund Monte Carlo programs: 

, No 
Jq{ x ) 



x(l - ± - ^2_)2 



where cq is a parameter dependent upon the heavy flavor and Nq is the corresponding 
normalization factor. 

The first postulation has been studied by several authors [11,12]. The universal distri- 
bution function of the QCD is interpreted as the structure function of the 6-quark in the 
parton model, which determines the distribution of the light-cone projection of the 6-quark 
momentum inside the B meson. In this and the next sections, we present the heavy quark 
expansion of the model in the viewpoint of the parton model. 

Being on a ground of an phenomenological model, the 6-quark mass mj is not an ingre- 
dient of the model and should be defined in the framework of the model in order to express 
the phenomenological model as l/rriQ expansion. We define the 6-quark mass as: 

m b = {x)m B , (3) 

following ref. [12]. Since we have the measured value of mg and well-defined function f(x), 
this relation explicitly defines the mass of the 6-quark. The nonperturbative informations 
of the parton model are encoded in the distribution function f(x). If we want to express the 
parton model with the heavy quark expansion, it is performed by expanding the distribution 
function in 1/rag. We write the expansion of the distribution function f(x) as 

f(x) dx = N ( 6(x B ) + — fi(x B ) + -Kf2(x B ) + ••• ) dx B (4) 
\ m b m( J 

where N is a normalization constant, which is expressed by the variable 



mB ( \ (<\ 

x B = ~^-\ x - x o) (5) 

chosen by Bigi et al. [9,11]. Here A is the effective mass of light degrees of freedom inside B 
meson and x = 1 — A/m B . The leading term is the Dirac delta function since we demand 
that the leading term of l/mj, expansion gives the free quark decay when the normalization 
constant is also written in the form of l/mj, expansion 

N=\l-\ / dx B fi(x B )-\ o / dx B f 2 (x B ) + ... ) 

\ m b J m 2 b J J 

Our definition of mj results in the well-known relation of meson and quark masses in 
the parton model. From the definition of the 6-quark mass, we obtain 



m B = m b + A / dx B x B f 1 (x B ) 

mj J 

-~2 (/ dx B x B f 2 (x B ) - J dx B x B f 1 (x B ) J dx B fi(x B 
+°(^) ■ ( 6 ) 



Provided we let 



A J dx B x B fi(x B ) = Am 2 B , (7) 



we obtain the mass relation 



— Am| 1 
m B = m b + A h 0{ — t 



m b m l 



III. INCLUSIVE SPECTRA AND THE DISTRIBUTION FUNCTION 

In this section, we explicitly relate the universal distribution function of the QCD to the 
parton structure function in order to help us to calculate the moments of the distribution 
function of 6-quark in the framework of the parton model. 

When we do not consider the perturbative QCD corrections, we can investigate the 
decay spectrum in cleaner way in the inclusive radiative decay process B — > ~{X S than in the 
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inclusive semileptonic decays. This process is not studied with the parton model framework 
yet. The relevant effective hamiltonian is given by 

Hen = -^V tb V t :c 7 (m b )-^—sa^(m b P R + m s P L )bF, v (9) 

in leading log approximation. Hereafter we neglect the s-quark mass. Since the partonic 
subprocess is 2-body decay, the photon energy spectrum is represented by a monochromatic 
line: 

dr<°) f 

_ = IV(£ 7 - -m b ) (fO) 

which will be the leading term of the inclusive decay spectrum of B — > ~{X S in 1 /m b expansion 
obviously. We write the decay width in the framework of the parton model 

E B dT(B 7 X S ) = J dx f(x) E b dT(b s 7 ) . (11) 

The decay rate of the partonic subprocess is directly obtained from amplitude, dT oc \Ai\ 2 = 
Gpe 2 \V tb V* s \ 2 \cr(rn b )\ 2 rnl / Att 4 since the amplitude is constant for 2-body decay. We have 
the photon energy spectrum 

dT 2 -r f( x+ ) (i + o(— )) (12) 



dE 1 ttib V m b ' 

where x + = 2E 1 /niB and r is defined in eq. (10). This spectrum is shown in Fig. 1. We 
note that the line shape is directly determined by the distribution function of 6-quark inside 
B meson as is observed in the Ref. [11,12]. 

In the framework of the HQET, Bigi et al. [11] derived the photon energy spectrum by 

dT 2 roF(x B )(l + 0( — )) (13) 



dE 1 A V m b 

where 

2 / 1 \ 

XB = A \ El ~ 2 m 7 ' ^ 14 ^ 
We find that the line shape is determined by the QCD distribution function F(xb) } which 
shows scaling behaviors of QCD without perturbative QCD corrections. However the genuine 
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function F(xb) cannot be determined completely by perturbative consideration. The HQET 
gives only some restrictions on this function, e.g. we can calculate moments of F(xb) in 
principle with the help of the HQET and OPE. Here we face the need of constructing realistic 
models with the help of informations from QCD. 

We find that the result of parton model also shows the scaling properties and that the 
distribution function f(x + ) is directly related the line shape function F(xb)- By some 
variables changes, we obtain the relation 

—/(*+) = i (—f(x + - x )) = If(xb) . (15) 

TTlB A \TTlB / A 

The function / is the one obtained by the variable change x + — > (x + — x ) + x from f(x + ). 
The functional form of F(xb) is obtained by the rescaling ( ) = (A/i71b)xb from 

f(x + — x ). Consequently the distribution function of QCD, F(xb) } is the very distribution 
function of 6-quark in our picture up to variables changes. 

We show that in the inclusive semileptonic decay process, the differential decay rate of 
the parton model is also directly related to that of the HQET after the appropriate variable 
changes. The scaling behavior in the inclusive semileptonic decay process is derived in the 
HQET, see ref. [11] 

dT G 2 F \V ub \ 2 2^^ .24(^0 - E^m^ - q 2 ) 

- r{x B ) (Id) 



dEidq 2 dq 1927T 3 A mj — q 

where 

m 2 h + q 2 - 2m b q 

X B = t7 ; 17 

Apart from kinematic factor, the triple differential decay rate is determined by one scaling 
variable xb- The scaling behaviors of QCD is violated by the perturbative and nonpertur- 
bative (higher twist) corrections. 

The parton model approach as is presented in eq. (1) gives the triple differential decay 
spectrum 

G^IK&I 2 qo — Ei 



dE t dq 2 dq 4tt 3 /q 2 + m 



(x + f(x + )(2Ei - m B x_) + (x + £_)) (18) 



where 

q ± Jq 2 + m 2 
x± = . 

m B 

Jin and Paschos [3] argued that the contribution of f(x_) is expected to be relatively small 
in the kinematic region where we are interested in and can be neglected by the kinematic 
analysis. This is in consequence of the fact the light-cone distribution function is sharply 
peaked around x ~ m^/m^. We let m q = hereafter. From the eq. (18), we obtain the 
coincident form with that of the HQET 

dT G 2 F \V ub \ 2 2 24( 9o - E,)(2m b E, - q 2 ) 

,n 7 2 j ~ mo 3 JV X + • l ly J 

dhjidq^dqQ iy27r J mg mj — g 

As the case of B — > ~{X S} the line shape function F(xb) is related to the distribution function 
f(x + ) of our model in the same manner after appropriate variable changes. We see that the 
parton model shows the same form of scaling behavior as that of QCD. It would also be 
violated by the perturbative QCD corrections and nonperturbative effects which are mainly 
related to the behaviors far from the light-cone. 

IV. MOMENTS OF THE DISTRIBUTION FUNCTION 

We have investigated that the first postulation of the parton model is consistent with 
that of the HQET. Now the validity of the model is attributed to the functional form of the 
distribution function f(x). Here, we study the second postulate, the validity of Peterson's 
fragmentation function as a distribution function of 6-quark. Using the HQET and the 
OPE, the moments of QCD distribution function 

dxs x b F(xb) , n = 0,l,..., (20) 

are calculated from the expectation values of local operators between B mesons: 
1 



2m B 



< B\Sbir lll ...Tr^b — traces|i? > = a n A n (v fll ...v fll — traces) (21) 
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where S denotes the symmetrization operator and tt^ = iD^ — m^v^. a = 1 from 6-number 
conservation, and a\ = up to (9(l/m 2 ) order indicates the lack of corrections. For 

the second moment of the spectrum, it is related to the average kinetic energy of 6-quark 
inside the B meson. We have 



With Peterson's function, the moments of the photon energy spectrum are numerically 
calculated. We have a = 1 from normalization. a\ is simply zero because we define the 



tub- In fact, our definition of 6-quark mass aims at this feature of QCD. The second moment 
is derived as 



where y = 2E 1 /niB- We obtain the numerical value a 2 — 0.757 with the parameter e = 0.004 
and the mass of B meson mg = 5.3 GeV, which well agrees with the result of the HQET and 
QCD sum rules, a 2 ~ 0.5 — 1. Thus one would conclude that the photon energy spectrum 
predicted from the parton model is in accord with the result of the HQET up to the order 
of 1/m 2 . This result explains the fact that the parton model present well-agreed spectrum 
with the experiment from ref . [2] . 

May we conclude that Peterson's form is appropriate as the genuine distribution function 
of 6-quark inside the B meson from this agreement? The answer is no. In our framework, 
the value of 6-quark mass is defined as mj ~ 4.5 GeV and the effective mass A ~ 0.8 GeV. 
These values show some differences from those predicted from the HQET, A ~ 0.4 — 0.6 
GeV. Moreover the average kinetic energy of 6-quark is calculated in the parton model 
/i 2 = 1.43 GeV 2 , which is too much larger than the value estimated from QCD sum rules 
[13] to be /i 2 ~ 0.6 GeV 2 . The former discrepancy means that Peterson's function gives 
too small mean value and the latter that it gives too large variance. These results make 
the conclusion that Peterson's function does not give the correct description of 6-quark 





6-quark mass is the average value of the distribution function f(x) times the meson mass 
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distrbution inside B meson comparing with that of real QCD. It gives allowable values of 
the second moment since above two effects are cancelled each other in calculation of a 2 . We 
will discuss the improvement of Peterson's function as a distribution function in section 5. 



We study Peterson's fragmentation function here. This function has been used as a 
description of fragmentation process of b- and c-quarks into hadrons containing one heavy 
flavour. The functional form is not purely ad hoc. It relates the energy transfer to the frag- 
mentation amplitudes. As the distribution function of 6-quark inside B meson, Peterson's 
function gives rather wrong values for some QCD calculations as shown in section 4. This 
function shows large asymmetric behavior with respect to the peak and it gives too small 
mean value which is too much deviated from the peak and too large width. We find that its 
functional form should be improved. 

The improvement is performed by adding a following physical argument, which makes 
better the functional form of the Peterson's function as an QCD distribution function. 
As stated before, Peterson's function describes the fragmentation process by the energy 
transfer of the process. In the vewpoint of the parton model, the principal features of the 
amplitude for a fast moving heavy quark Q fragmentation into a hadron H and a light 
quark q are determined by the value of the energy transfer AE = Eh + E q — Eq such that 
amplitude(Q — > H + q) oc AE -1 . On the other hand, we are considering the decay process 
in this paper. If we accept the assumption that the amplitude is determined by the energy 
transfer, the energy transfer should be AE = Eq + E q — Eh in decay preocess. This is 
the same as the fragmentation process when we let mg = mj as Peterson et al. assumed 
[10]. We know that actually mg > mj and we introduce a new parameter ctg to express the 
mass difference such that m B = m\{l + olq). With this alteration, the energies about the 
transverse particle masses are given by 



V. IMPROVEMENT OF PETERSON'S FUNCTION 
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-« + p 2 ) 1/2 



1 

oc 1 + a Q 



1 - z 



e 



(24) 



where z is a ratio of 6-quark momentum to B meson momentum. Then we obtain the 
function 



We show the functional form of this function in Fig. 2 compared with the original Peterson's 
form with the parameter cq = 0.004 and aQ = 0.085. As we expected, this new function 
gives larger mean value and smaller variance. 

Using this improved function with the parameter cq = 0.004, aQ = 0.085, we obtain 
the numerical value for a few moments. Two moments a and a\ are not changed. The 
second moment a 2 = 0.673, which gives the proper line shape. We also obtain the values of 
the effective mass A = 0.537 and fj 2 v = 0.583, which are consistent with those obtained by 
the HQET and QCD sum rules. Up to this order, the parton model with our improvement 
gives the consistent results with those of the genuine QCD distribution function estimated 
from QCD sum rules. We present the moments and the values of A and fj 2 v with varying 



We study the parton model with the heavy quark expansion through the inclusive radia- 
tive decay and the inclusive semileptonic decay of B mesons. We properly define the mass of 
6-quark in the framework of this model and write the expansion of the distrbution function 
of 6-quark in l/m^. By relating the 6-quark mass to the meson mass with this heavy quark 
expansion, we see that our definition of the 6-quark mass satisfies the features presented by 



In the parton model, the distribution function of 6-quark, which describe the motion of 
the quark inside the B meson, plays a major role. It has been shown that the distribution 




(25) 



a Q =0.02-0.1 in Table 1. 



VI. CONCLUDING REMARKS 



QCD. 



11 



function f(x) is directly related to the QCD distribution function determining the line shape 
of the decay spectra. Here we use Peterson's fragmentation function as distribution function 
of 6-quark following Ref. [1,2]. However the original Peterson's function gave rather wrong 
values for the average kinetic energy of 6-quark inside B meson and the effective mass of 
the light degrees of freedom while it gave allowable numerical values for the second moment 
a 2 . Thus improvement of the function is required. We observe that Peterson's function 
assumes that mj = mg. By parametrizing the mass difference of 6-quark and B meson, we 
rederived the distribution function for decay process and obtained the refined form eq. (24). 
We showed that this function satisfies the restriction of QCD and conclude that the model 
is successfully improved. 

In Fig. 3, we compare the decay spectra calculated by Peterson's function with those by 
our refined function. We find that the lepton energy spectrum of B — > evX u decay does not 
show mush change in spite of the drastic change of the distribution function. This resulted 
in the fact that Jin et al. could obtain well agreed spectrum with experiments by using 
Peterson's function as the distribution function in ref. [2]. 
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